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An exact solution to the dynamic equations for a massive boson traveling in a pp-wave gravitational background 
under the influence of the force induced by curvature, is presented. We focus on the effect of anomalous polarization- 
curvature interaction and consider models in which the coupling constant of such an interaction is treated to be either 
a deterministic quantity or a random variable. 



1 Introduction 

The motion of a point particle with intrinsic structure has 
become a subject of detailed studies about a century ago. 
M. Abraham £Q was the first to consider in detail the mo- 
tion of an electron as a particle possessing a supplemen- 
tary vector degree of freedom. In the twenties of the past 
century L.H. Thomas [2], J. Frenkel 3 and I. Tamm @| 
studied different aspects of the classical spinning parti- 
cle motion. The equations of motion for particles with 
arbitrary multipole moments in the framework of Gen- 
eral Relativity were obtained by M. Mathisson and 
then revived by A. Papapetrou B] and W.G. Dixon [7]. 
Studying the effect of spin precession in an external elec- 
tromagnetic field in the framework of special relativity, 
V. Bargmann, L. Michel and V.L. Telegdi introduced 
the term "anomalous magnetic moment of an electron" . 
In the absence of an anomalous moment, the Bargmann- 
Michel- Telegdi (BMT) equations coincide with Bloch's 
equations 9]. In the last three decades, the interest in 
studying the dynamic equations for particles with intrin- 
sic degrees of freedom, being derived in the framework 
of the Lagrangian or Hamiltonian formalism, has grown 
considerably [10-33]. Thereupon, the BMT-model has 
gained the status of a test model for verification of new 
models of dynamics of a point particle with spin or polar- 
ization in external fields of various kinds. The history of 
studying the dynamics of spinning particles is presented 
in the book |34j : alternative approaches to derivation of 
motion equations for charge and spin, the problem of non- 
collincarity of momentum and velocity four-vectors and 
the origin and interpretation of the Thomas precession 
are examined in detail there. 

The most important thing while analyzing the dy- 
namic equations for particles with intrinsic structure in 
the gravitational field is the studying of a spin-curvature 
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(polarization-curvature) interaction. Per se, the Mathis- 
son-Papapetrou-Dixon equations are a non-minimal gen- 
eralization of the classical dynamic equations because 
they take into consideration the Riemann curvature ten- 
sor and its covariant derivatives. Formally, the introduc- 
tion of curvature-dependent forces (or "tidal" forces by 
the conventional terminology) breaks the free fall uni- 
versality principle and thus contradicts Einstein's equiv- 
alence principle. Thereupon, the spin (polarization) - 
curvature interaction is now considered to be a theoreti- 
cal prediction whose experimental verification may have 
an interesting application to the problem of testing the 
metric theories of gravity. Following LB. Khriplovich 
E01 12SI I2H EDI and M. Bander and K. Yee |2U, using 
an analogy with the anomalous interaction of a spin- 
ning particle and the electromagnetic field, introduced 
in the BMT theory, it seems reasonable to consider 
also the anomalous interaction of spin (polarization) 
with the gravitational field. Such an analysis involves 
a new parameter Q as a constant of anomalous spin 
(polarization)-curvature interaction, by analogy with the 
phenomenological parameter g — 2 , introduced in the 
BMT theory [H] as a constant of anomalous interaction of 
spinning particle with the electromagnetic field. Under a 
generalization of such kind, spin (polarization) dynamics 
is predetermined by a Lorentz-type force, one addend 
of which is proportional to the Riemann curvature ten- 
sor pi) . 

In this paper we suggest an exactly integrable model 
of the dynamics of a relativistic particle with an arbitrary 
directed polarization vector in the gravitational radiation 
field. The model is based on the known dynamic equa- 
tions for particles with supplementary degrees of freedom 
and takes into consideration an anomalous interaction of 
polarization with curvature. The paper develops our re- 
cent works 1201 OH EH] ) m which, in particular, the 
effect of parametric oscillations of Faraday rotation of 
the polarization vector was predicted. The results ob- 
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tained in this paper demonstrate a possibility of new type 
for the precession of the polarization vector, indicated as 
"hyperbolic rotation" . We consider two types of mod- 
els: first, a model with a deterministic coupling constant 
of anomalous interaction between the space-time curva- 
ture and particle polarization, and, second, the model 
with a random coupling constant. One of the purposes of 
this paper is to compare the behaviour of particle in two 
different cases: first, when the anomalous polarization- 
curvature interaction is absent (the coupling constant is 
zero identically), second, when the random coupling con- 
stant is zero on the average. 

The paper is organized as follows. In Sec. 2, evolu- 
tion equations for the velocity four-vector and the spin 
(polarization) four-vector are established and compared 
with the known dynamic equations for particles with an 
intrinsic structure. In Sec. 3, exact solutions to the re- 
sulting evolution equations are found for the case when 
a boson travels in a pp-wave gravitational background 
under the influence of a deterministic force induced by 
curvature. In Sec. 4, we discuss the specific features of 
boson motion in the case when the coupling constant of 
the polarization-curvature interaction is considered as a 
random variable with zero mean value. The last section 
contains conclusions. 



2 Evolution equations 
2.1 General model 

It is well-known that, for a point particle with supplemen- 
tary tensor degrees of freedom, the generalized momen- 
tum four-vector P % is not parallel to the velocity four- 
vector U % = (see, e.g., 151) 1 



P l = mcU 1 + q 1 , q r U % = . 



(1) 



The part of the momentum q % , which is orthogonal to 
the velocity four-vector, may include the angular mo- 
mentum tensor S tJ , the quadrupole momentum tensor 
jtmns ag we jj ag tensors describing higher multipole 
moments [51 128| 



(2) 



The dot denotes the derivative D/Dt = U l Wi, where 
Vj is the covariant derivative. The Riemann curvature 
tensor R nm sj and its covariant derivatives happen to be 
included into the formulae for the generalized momen- 
tum as necessary structure elements. As usual, the scalar 
product UiU 1 = const is set equal to unity identically, 



UJT = 1 . 



(3) 



The scalar product P % Pi = m 2 c 2 + q l qi does not remain 
constant along the particle word-line, or, in other words, 
it is not an integral of motion. In our opinion, it is 
then more convenient to construct the dynamic model 



for particles with intrinsic structure using the velocity 
four-vector then. Starting from the relation 



U i U i = ^{U i U i ) = Q, 



(4) 



one can find the generic form of an evolution equation for 
the velocity four-vector: 



U l =f\ fu,=o. 

Using the identity: 
fi = Ui(f k U k ) + ^f k , 



A? = 5? - U l U k 



(5) 



(6) 



and the orthogonality condition f k U k = 0, one can rep- 
resent the effective force four- vector /' in the form 



f = A k V k {x, U, S) = (ViU k - U z V k )U k 



(7) 



Here V k (x, U, S) is an arbitrary vector function depend- 
ing on the coordinates, the velocity four-vector and spin 
(polarization) variables. Introducing the new antisym- 
metric tensor (simple bivector) 



Cl ik (x,U,S) = V i U k -U i V k , 



(8) 



one can represent the dynamic equation ||S} in two equiv- 
alent forms 



W = Aj V k (x,U,S) 



U% — Qik{x, U, S)U k . (9) 



To characterize the state of a particle which possesses a 
vector degree of freedom (spin or polarization) , one uses 
the well-known decomposition of the anti-symmetric ten- 
sor of total moment 



Sik = LiU k — L k Ui — eikimS U 



l t jin 



(10) 



Here L 1 is the four-vector of orbital moment and S l is 
the spin (polarization) four- vector, defined as follows: 



Li — Si k U k 



<? = 1 



where e ik i r , 

e%klra — * 



is the Levi-Civita tensor 

—gEiklrn , E0123 = — 1 , 



(11) 



(12) 



and the asterisk in Ijllfl introduces the corresponding du- 
ality symbol. The four-vectors Li and Si are orthogonal 
to the velocity four- vector automatically by the definition 

CD 



L l Ui = 0, S k U k = 0. 



(13) 



In this sense the so-called Tamm condition S k U k = 
is not an additional requirement. One can find an ob- 
vious analogy between the decomposition l|ll(l and the 
decomposition of the Maxwell tensor F% k in electrody- 
namics [55] 



Fik = EiU k — E k Ui — eikimB U 



Irjin 



(14) 



Following this analogy, one can say, that S l plays the 
same role as the four-vector of magnetic induction B l 
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(and is in fact a pseudo- vector) , while L 1 plays the role 
of electric field four- vector E % (and is a true vector) . The 
four-vector S l is considered to be a vector parameter of 
the intrinsic state of a particle. The scalar square of this 
four-vector is assumed to be constant along the particle 
world- line, i.e., 



S k S k = -E 2 = const ^ . 
This requirement results in 

s k s k = , 



(15) 



(16) 



and, following the previous discussion concerning the U l 
decomposition, one can state that the corresponding evo- 
lution equation has two equivalent forms 



S l = 5? 



S t S k 
SiS 1 



W k (x,U,S) < ► 

S l = io lk (x,U,S)S k , 



where 



WjS k - s t w k 
s l s. 



(17) 



(18) 



From the orthogonality condition SiU 1 — it follows that 

U l S l + S'U l = 0, (19) 

consequently, the four-vectors Vi and W% as well as the 
tensors £li k and Wi k are connected by the requirements 



u l w t + sty = o , (io lk - n lk ) u i s k 

This means, in particular, that 

uj lk s k = n ik s k + n lk {x, u, s)u k , 



0. 



with 



-n 



kt • 



VL ik U i S k =0. 



(20) 
(21) 
(22) 



Thus, we can represent the evolution equations of a parti- 
cle with spin or polarization in the following generic form: 



Ui = n lk u k 



Si 



fli k S k 



(23) 



2.2 Comparison with well-known models 
2.2.1 The Bargmann-Michel-Telegdi model 

To recover the BMT results [H] 



Si = 



9 jP Qk 



U 



9 



kl >- 



F lk U k 



-l)F M S k U l 



(24) 



where g is a gyro-magnetic ratio, one has to suppose that 
Qik is proportional to the Maxwell tensor Fi k : 



Qik = 



ik 



and Oife has to be equal to 



Qik = 



(l-l)(U k F a -U iFkl )S> 



(25) 



(26) 



2.2.2 The Bander - Yee model 

Bander and Yee obtained the following spin evolution 
equation: 



DS ab 

Dt 

( eg 



{UbSac — U a Sbc) 
K 



DU C 



■F 



cd 



R cdef S, 



Dt 

'ef 



\ Amc 2 8mc 

x{SacA bd + S bd A ac ). (27) 

The dynamic equations derived in [52] have terms which 
are second-order in the spin, terms with derivatives of 
the Maxwell tensor and a term with the second order 
derivative of U % . In the approximation where we neglect 
all these terms, we have the following equation |24|: 

1 



Ui = 



tF 



ik 



2mc 



tj cat 
J^ikab^ 



u« 



(28) 



Let us consider a model with vanishing orbital moment 
(Li = 0). Then, in terms of the spin (polarization) four- 
vector 5* a , the dynamic equations l|28|) have the form (|23|l 
if 



where 



; Fa 



1 



R*ikabS a U 



D* p cd 

K-ikab — 2 ik e cdab 



(29) 



(30) 



is the right-dual Riemann tensor. This quantity is a 
pseudo-tensor, but the product R* kab S a happens to be 
a true tensor. Eqs. (|27Jl can be represented in the form 
(J23) when 



wit i 

^ik = 



u k riu — f/jOfe 



s l . 



:F 



ik 



(§-') 



hi — 1 



mc 



RlkabS^ 



(31) 
(32) 



2.3 Equations of boson dynamics 



Let us consider the dynamics of vector bosons possessing 
zero charge e = and polarization four- vector S l . The 
dynamic equation for such a particle can be immediately 
obtained from 



U 



1 



~R*iki m u k s l u m 



(33) 



The equation of polarization evolution, obtained from 
(PI) with (23, Ell) and (EH, yields 

Si = ~ [KR*k lm s k s l u m 

+ ( K -l)UiR* klmn S k U l S m U n ]. (34) 

Eqs. ESI and EH can be also obtained from the BMT 
equations by a formal substitution discussed by Khriplo- 
vich in |30j : 



-Fik 
c 



p Qa b p* carrb 



(35) 
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if one replaces g/2 by n. One sees that the case k = 1 
corresponds to the BMT model without an anomalous 
magnetic moment. We introduce a new dimensionless pa- 
rameter Q by the definition Q = k— 1 and suppose that, 
from a phenomenological point of view, the Q parameter 
is a coupling constant describing an anomalous interac- 
tion of polarization with curvature. Finally, one obtains 
the following equation of evolution for the polarization 
four- vector: 

S i = --±-[(Q + l)-R* klm S k S l U m 



+Q ■ U t Rt lmn S k U l S m U n ] 



(36) 



3 Boson dynamics in a pp-wave gravitational 
field. Deterministic model 

3.1 Gravitational wave background 

The well-known pp-wave solution of the vacuum Einstein 
equations can be represented by the metric |31| 



+ e 



where 



-2/3 



1 



dr 2 = 2dudv - L 2 { [e 20 (dx 2 ) 2 
(dx 3 ) 2 ] cosh 27 + 2 sinh 2^dx 2 dx 3 } 



(ct-x 1 ), 



1 



(ct 



(37) 



(38) 



are the retarded and advanced time, respectively. The 
functions (5{u) and ^(u) are considered to be arbitrary 
functions of the retarded time, satisfying the following 
conditions on the wave front plane u — : 

0(0) = 7 (0) =0, /3'(0)= 7 '(0) = 0. (39) 
The background factor L{u) satisfies the equation 

L" + L (f3' 2 cosh 2 2 7 + 7 ' 2 ) = (40) 
and the initial data 

L(0) = 1, 1/(0) =0. (41) 

The gravitational wave (GW) is indicated as a GW with 
the first polarization when (3{u) 7^ 0, 7(1*) = 0, and 
as a GW with the second polarization when j3(u) = 0, 
j(u) ^ 0. For the pp-wave solution lpT7l) . the right- 
dual Riemann tensors R*j Mm (|3*Uj) and the left-dual one 



x jklm — 2 c jfc 

nonvanishing components 



Cau 9 Rpqim coincide and have the following 



rL 2u2u 



L 2 ^ 7'- -/?'sinh4 7 



R 



3u3u 



L 2 e 



-20 



2' 

7' + -0 sinh 47 



R*2u3u=R*3u2u= [£ 2 /3'cOSh 2 2 7 ]\ 



(42) 



Here and in what follows, Latin indices run over four 
values: u, v, 2, 3, while Greek indices run over two values: 
2,3. 



3.2 Exact solutions to the evolution equations 

3.2.1 Key subsystem of the evolution equations 

The GW metric l|37|) admits two spacelike Killing vectors 
and one null covariantly constant Killing vector |31| 



e=si, e=si, e=K, 

(2) (3) (v) 

9lk ee = o, v fc f = o. 

(v) (v) (v) 



(43) 
(44) 



The projection of the dynamic equation 1)33(1 onto the 
direction given by £ l yields 



(v) 

— — > U v = const = C v 

dr 



(45) 



It follows from (|45|) that the parameter r is linearly con- 
nected with the u coordinate since 



^^u u = u v = c v 

dr 



(46) 



When C v (this condition does not hold only for a 
massless particle co-moving the GW) , one obtains 



To 



C11 



(47) 



Since U t Ui = 1, the U u component of the velocity four- 
vector can be found immediately: 



U u = ^ [1 - g aa U a U a ] 



(48) 



Similarly, from S l Ui — 0, one obtains an expression for 
the S u component of the polarization four-vector: 



S u = --^[S v U u + g a °U a S a } 



(49) 



Thus the U u and S u quantities are found as soon as U2 , 
U3 , S2 , S3 and S v become known. For the remaining five 
components of the velocity and polarization four-vectors, 
one has five equations: 



dUg __ J_ » 

au mc 

dS n 



{C V S° - S V U°) , 



C v — : — (C v S a — S v U a ) 

du 

9aa _ Q 1 r. r>* 

2 mc 



— U a R* puau (C v Sr ~ S V U») 
mc H 



dS v 
du 



Q 



T? ■ ■- 



(C v S a - S v U a ) 
x (CvS° - S V U°) 
with the following initial conditions: 

U a (0) =C a , S a (0) = E a , S v (0) = E v 



(50) 
(51) 
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After introducing the new variables 
X a = C v S a S v U a 



(52) 



a closed subset of three equations emerges from the 
set JSD) 



dX a 
du 

dS v 
du 



where 



a°(u) - —S v R*° uau (u) 
mc 



Q 



Mr, 



t (u)X<*X° 



<{u) = -g ap (u)g' pa (u). 



(53) 
(54) 

(55) 



The prime denotes a derivative with respect to the re- 
tarded time u. For the transversal components of the 
velocity four-vector U a we obtain two decoupled equa- 
tions 



dU Q 1 

— rt n 



du mc 

which give the following formal solutions for U a : 



U a {u) = C a 



1 

mc 



I! 



t (u)X a (u) du. 



(56) 



(57) 



The unknowns S a (u) can be found immediately from the 
relation J52J): 



S a (u) = — [X a (u) + S v U a (u)} 



(58) 



Thus all the solutions of Eqs. I|50|l will be represented 
by quadratures as soon as solutions to the system 
(|54|l are known. We will designate the latter as the key 
subsystem. 

3.3 Exact solutions to the key subsystem for the 
case Q = 

One of the purposes of the paper is to compare the be- 
haviour of a vector boson in the GW background in two 
cases. The first one takes place when the coupling con- 
stant Q is zero identically, i.e., when there is no anoma- 
lous interaction of polarization with curvature. The sec- 
ond case is connected with the suggestion that Q effec- 
tively depends on the environment of the boson and is 
a random variable with zero mean value. Thus we first 
consider the evolution equations with Q = . From (|54|l 
it directly follows that 



S v (u) = const = S v (0) 
and Eqs. H53[) result in 
X' a =a a a {u)X (T . 



E„ 



(59) 
(60) 



Such a system have been solved in |40| and discussed in 
[35-38,41,42], and here we use the result: 

X 2 {u) 
X 3 (u) 

cosh 7 sinh.7 
J \ sinli7 cosh 7 

cosV' —simp \ ( C V E 2 — E V C 2 
simp cosip / I C V E 3 — E V C 3 



(61) 



where 

u 

ip(u) = J (3' (u) sinh 2 j (u) du. 





(62) 



The precession of the polarization four-vector, described 
by Eqs. ((SHI, EH, (E2J> was discussed in detail in 

Ref. (23 . 

3.4 Exact solutions to the key subsystem for the 
case Q 7^ 

To gain analytic progress in the case of an anomalous 
interaction of polarization with curvature, let us consider 
the GW to have the first polarization (7 = 0). The right- 
dual curvature tensor R* klmn (|4*2"|l has only two nonzero 
components: 



and thus the equations for X a yield: 



(63) 



X 2 . (65) 



X' 



V 



L J mc " L 2 

The substitution 

X 2 = Le p Y 2 , X 3 = Le-PY 3 



(66) 



results in the following equations for the new variables Y 2 
and Y 3 : 



y{ = n(u)Y 3 , y 3 ' = n( u )Y 2 , 

Q [L 2 (u)P'(u)]> 

Sl{u) = S v {u) —— 

mc L z (u) 



(67) 
(68) 



Solutions to this system can be represented in the ele- 
mentary functions: 

Y 2 {u) = Y 2 {0) cosh*(u) + Y 3 (0)sinh*(u), 

Y 3 (u) =F 3 (0)cosh*(ii)+y 2 (0)sinh*(w), (69) 



where 

* = / n(u)u, *(o) = 0, 

Jo 



(70) 
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y 2 (o) = C V E 2 - E V C 2 , 
Y 3 (0)=C v E 3 -E v C 3 , 



(71) 



As for the unknown S v (u) , using the expressions for 
quadratic integrals of motion U l Ui=\, S l Ui—0, S % Si——E^ 
as well as the definition of X a (u) 1)52(1 . and the expres- 
sions for Y a (u) H69fl . one can show that the following 
relations take place: 

S 2 v (u) = C 2 v El - [Y 2 2 (0) + Y 2 (0)] cosh2*(u) 

-2y 2 (0)y 3 (0)sinh2*(u) = B 2 (^) , (72) 

E 2 =C 2 E 2 -Y 2 (Q)-Y 2 (0). (73) 

The relations (EEJ) , (EOJ > make it possible to find *(w) 
by quadratures: 



*(«) 



± 



o 



mc 



where 



r \u)= / R 3 u3u (u)du 



(L 2 f3>) 
L 2 



(74) 



(u)du . (75) 



The function r(u) has the dimension of distance; it may 
be considered as some effective curvature radius. 

Generally, the solution for \&(it) can be represented in 
terms of elliptic functions. Below we consider some ap- 
proximations to such solutions. When the function ^(u) 
is known, the transversal components of the velocity four- 
vector and the polarization four-vector take the form 

U 2 {u) = C 2 + ^^- J S [ [L 2 f3'}'^-cosh^(u)du + 

U 

J /[L^'^sinh*! (u)du, (76) 



o 

^2(0) 

mc 



U 3 (u) = C 3 + 



1^(0) 



mc 



[L 2 (3'}' cosh* \ (u)du + 



o 



-0 



\L 2 (3'}' — sinh * } (u)du , (77 1 
mc J l L ' 

o 



S 2 (u) 



^LMe^W cosh*(?i) 



+YM L{u y(u) M(u) + Sv[u) ^M , (78) 

S 3 (u) = ^I(ii)e-"W cosh*(u) 
+ YM L{u)e -P(u) sinh v(u) + S v (u)^± . (79) 



Eqs. ((76(1 - 1(79(1 display a new type of precession of the 
polarization vector, which can be designated as hyper- 
bolic rotation of polarization under the influence of an 
anomalous curvature force. The above quantities depend 
non-linearly on the Q parameter via the function \I/(it) 
given by Eqs. and J72JI. When Q = 0, = 0, 

and the integrals in 1(761) and 177fl take an explicit form: 

02|Q=o(«) = C 2 + —{C V E 3 - E v C 3 )[L{u)S u) \ , (80) 



U 3lQ=0 (u) =C 3 (C V E 2 - E v C 2 )[L(u)e-^}' , (81) 



,/3(«) 



&-^g][L(«)e™r, (82) 



S 3 | Q=0 ( M ) = + (# 3 - ) L(u)e^> ' 



■^(e 2 -e v ^ 

mc \ C„ 



[L{u)e 



(83) 



Thus Eqs.gHJ), (EH, JZnj-JEII) describe the 

evolution of a relativistic polarized boson in a pp-wave 
gravitational background in the framework of a model 
with the deterministic parameter Q (coupling constant of 
the polarization-curvature interaction). Eqs. ((80(1 - 1(83(1 . 
describing the particular case Q = , will be used directly 
in the next section. 

4 Random polarization-curvature interaction 

We suppose that the anomalous interaction of polariza- 
tion with curvature has a random nature and can be at- 
tributed to the particle interaction with its environment 
mediated by the space-time curvature. Mathematically, 
this means that the parameter Q can be treated as a ran- 
dom variable. Basing on such an ansatz, one can calculate 
the mean values of the velocity and polarization four- 
vectors and compare the results with those for Q = 0, 
given by Eqs. I(80|) - ((83(l . The distribution function for 
the Q variable is supposed to be Gaussian: 



f(Q) = 



i 



ex P 



Q_ 2 

v 2 



(84) 



where I? is a dispersion parameter. 

4.1 First special solution: Y 2 (0) = Y 3 (0) = 

To obtain a starting exact result, let us suppose that both 
initial values 1^(0) and 1^(0) are equal to zero. In this 
case it follows from (JT2J that S v (u) = E v = and 
*5f(u) can be found explicitly: 



*( u ) = -Q^r-\u) 
mc 



(85) 
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On the other hand, the ^(u) function happens to be a 
hidden one, since in this degenerate case U a , S v and S a 
remain constant: 



X a (u)=0, U a {u) = C a , 

E 

S(x{u) — —C a , S v (u) = E v = C V E . 
The quantities S u (u) and U u (u) 
S u (u)=-^[l + g^(u)C a C fi ], 



(86) 



(87) 



do not contain ^>(u) and consequently do not depend on 
Q . This special type of particle motion is characterized 
by the following feature: the projection of the polariza- 
tion four-vector onto the GW front plane is parallel to 
the projection of the velocity four-vector not only at the 
initial moment, but during the whole period of evolution, 
too. 

4.2 Second special solution: E v = 0, 

F 2 (o)-y 3 (o) = o 

According to l|53|) . Eqs. I|61|l . (|62|l describe an exact so- 
lution to the evolution equations also in the case when 
Q / 0, but S v (u) = 0. It is possible, when, first, 
S v {0) = E V = 0, second, X 2 {0) ■ X 3 {0) = C 2 U E 2 E 3 = 0, 
as can be seen from and In particular, it is 

possible, when the initial polarization four- vector S l (Q) 
has a vanishing projection onto the null Killing vector £ l 

(«) 

as well as onto one of the spacelike Killing vectors f 1 or 

(2) 

. The polarized boson with such characteristics, travel- 

(3) 

ing in the GW background, does not feel the influence of 
the curvature-induced force. Both in first and the second 
special cases, exact results of averaging over Q coincide 
with those for Q = . 

4.3 Third special solution: Y 2 (0) = -Y 3 (0) ± 

Eqs. lEU) show that if F 2 (0) = -F 3 (0) = Y(0), then at 
an arbitrary retarded time instant 

Y 2 (u) = -Y 3 (u)=Y(0) e~*W 

Then Eq. (|74l takes the form 

± I . 

>C 2 E 2 -2Y 2 (0)e- 21 - 



— r'^u). (89) 
mc 



Integration in (|89|) yields 



-*(«) 



C v Eq 
V2Y(Q) 



x cosh 



Ct,i?o _ QC V E _ x 

Arcosh— = =p r («) 

V2F(0) mc 



(90) 



At the initial instant u = 0, when r x (0) = 0, one ob- 
tains \l/(0) = 0. Then at the instant u* , when 



C V E QC V E Q _ x 

cosh r [u ) . 



(91) 



V2Y(0) 

the function c~* obviously reaches its maximum value 

_-*(«*) _ 



\/2F(0) 



(92) 



The instant u* is a stopping point for the hyperbolic ro- 
tation since the quantity S v (u*) vanishes (see, JZ2J)), and 
we obtain Q(u*) = and ff |u=u » = 0. Using in 
Eqs. I|76|) - (|79l) . we obtain an example of an explicit rep- 
resentation of an exact solution by elementary functions. 



4.4 Fourth special solution: Y 2 (0) = F 3 (0) / 

Exact results for this case can be obtained from the re- 
sults of the Sec. 4.3. by the formal substitution 'f — > — W . 

4.5 Statistical averaging in the general case 

Eqs. |(7Hl and ((74^1 show that, generally, the dependence 
of \& on the parameter Q can be represented in terms 
of elliptic function. For special initial data (see Sec. 4.1- 
4.4) one can obtain the results in elementary functions. 
Taking into account the relation l|72(l . one can state, that 
the growth of *&(u) is restricted by the requirement S 2 > 
0. When S 2 — 0, the function ^/ reaches its extreme 
values 



exp{2 *^ 2} = (r 2 (o) + y 3 (o)) 2 • (93) 

Both extreme values V^i and ^> 2 correspond to the rela- 
tions 



n(u*) = 0, 



du 



0. 



(94) 



When F 2 (0) -> -Y 3 (0), the expression §53$ gives JUJ, if 
we choose the minus sign in the numerator; the plus sign 
gives the solution for the case Y 2 (0) — > I3 (0) . Thus the 
function ^(u) remains finite. 

To perform averaging over Q , we consider the leading- 
order terms in the expressions for U a and S a . The di- 
mensionless parameter (E v / mc)^ 1 (u) is considered to 
be small. Since E v in quasiclassical theory is supposed 
to be proportional to the Planck constant, the inequality 
(E v /mc)r~ l {u) « 1 does not require, in general, that 
the GW is weak. When Y 2 (0) and Y 3 (0) are not zero si- 
multaneously, the leading-order term in the expansion of 
the function ^>(u) has exactly the same form as in Q85fl. 
Averaging with the Gaussian function Ij84(l yields 



= V 



E 2 v Y 3 (0) 
4m 3 c 3 



(U 2 )q - U 2{ q=o 

:(L 2 p>y 



L 



(u)dii . 



(95) 



Anomalous Polarization- Curvature Interaction in a Gravitational- Wave Field 



105 



V 2 E 2 V Y 2 (0) 
4m 3 c 3 



(S v }q — S v Q =0 



(U3)Q - C/3,Q=0 

-b(L 2 P'Y - 2 \r~^~ 
, p \ z y ( u )du, 



V 2 



(96) 



, o • T 2 (u) 

2E v m 2 c 2 

x{2K, 2 (0)F 3 2 (0) + E 2 v [Y 2 2 (0) + F 3 2 (0)]} , 



(97) 



(S 2 )q - S 2i q 

, V 2 YMEl 
4m 3 c 3 C„ 



V 2 Y 2 (0)E 2 _ 2 



Am 2 c 2 C v 

(L 2 p)'^-r-- }l<i)d„ 





P 2 {2r 2 2 (0)r 3 2 (0) + e 2 [y 2 (o) + r 3 2 (o)]} 



2m 2 c 2 E v C v 
mc 



r- 2 (u), (98) 



(S 3 )Q - S 3 ,Q=o = TZZTjTF, Le r W 



Am 2 c 2 C v 



l p 2 y 2 (o)-Eg 

D 2 F 2 (0)y 3 2 (0)^ _ x 



(L 2 /3')' 



(u)du 



7 3 C^Cy 



r- x {u) / <^ (L 2 f3'Y-—r- 1 \ (u)du 



L 



V 2 {2Y 2 (0)Y 2 (0) + E 2 [Y 2 (0) + Y 2 (0)]} 



C 3 

We have used the notation 
(A) Q = / dQf(Q)A(Q) . 



2m 2 c 2 S t ,C„ 



r~ 2 {u) . (99) 



(100) 



Thus the effect of random polarization-curvature force 
can lead to essential changes in the particle state when 
the dispersion parameter T> is large enough. 

5 Conclusions 

1. We have formulated the master equations govern- 
ing the dynamics of a massive relativistic boson with 
an arbitrary directed polarization four-vector, traveling 
in the gravity field under the influence of an anoma- 
lous polarization-curvature interaction (see Eqs. and 
(|36l) h The master equations constitute an analogue of the 
Bargmann-Michel-Telegdi model |Sj , describing the evolu- 
tion of an electron with an anomalous magnetic moment. 



They use the direct analogy between the electromagnetic 
and curvature interactions proposed by Khriplovich 30 , 
and represent a submodel of the Bander and Yee general 
model |2l) . 

2. The master equations have been solved exactly 
for the model of boson dynamics in a pp-wave gravita- 
tional background. We have distinguished three exactly 
solvable submodels. The first submodel describes the de- 
terministic anomalous polarization-curvature interaction 
(see, 63, (©, 6HJ), O, f7B | -f75 | ). The second one 
represents its particular case with a vanishing coupling 
constant of the anomalous polarization-curvature inter- 
action (see, l|80 |l -(|X3 |l 'l. The third submodel introduces a 
random anomalous polarization-curvature interaction. 

3. The exact solutions of the master equations ob- 
tained here demonstrate a specific type of the polariza- 
tion four-vector behaviour, designated as "hyperbolic ro- 
tation" induced by anomalous interaction of polarization 
with the space-time curvature. "Hyperbolic rotation" of 
the polarization four-vector affects the boson dynamics; 
the velocity four-vector depends non-linearly on the ini- 
tial data describing the direction of the polarization four- 
vector. 

4. The master equations happen to be nonlinear with 
respect to the velocity four-vector U % as well as the po- 
larization four- vector S l . This feature makes essentially 
different the two models: the first one, where the de- 
terministic anomalous polarization-curvature interaction 
vanishes, and the second one, where the random coupling 
constant of the anomalous interaction vanishes on the av- 
erage. The corresponding distinction is demonstrated by 
Eqs. JHSJ-S- 
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